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I«  The  Aim  of  This  Research 
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A  great  deal  of  the  existent  knowledge  of  the 
hydrogen  bond  in  the  solid  state  has  been  derived  from  X-ray 
crystallographic  studies.  The  solution  of  the  X-ray  crys¬ 
tal  structure  of  sodium  bromide  dihydrate  will  add  to  this 
basic  information  concerning  the  hydrogen  bond. 

It  is  of  interest,  at  this  point,  to  briefly  sum¬ 
marize  the  work  which  eventually  led  to  this  thesis  in¬ 
vestigation. 

It  has  been  reported  (Mutter,  Mecke  and  L^Ittke, 
1959)  that  a  hitherto  unknown  form  of  ice  resulted  when  water 
vapor  was  slowly  condensed  onto  a  sodium  chloride  plate  kept 
at  a  low  temperature.  The  infrared  spectrum  of  this  sub¬ 
stance  showed  sharp  lines,  in  marked  contrast  to  the  spectra 
of  cubic  and  of  hexagonal  ice.  A  repetition  and  extension 
of  this  work  at  Princeton  University  (Hornig  and  Schiffer, 
1961)  indicated  that  the  proposed  new  form  of  ice  is,  in 
reality,  the  dihydrate  of  sodium  chloride.  Comparison  of 
the  infrared  spectrum  of  sodium  chloride  dihydrate  with  that 
of  sodium  bromide  dihydrate,  produced  by  condensation  of 
water  vapor  upon  a  sodium  bromide  plate,  suggested  that  the 
two  structures  are  very  similar,  Schiffer  and  Hornig  also 
studied  these  salt  hydrates,  again  by  Infrared  methods, 
after  having  Introduced  small  replacements  of  HOO  for  the 
H2O  or  D2O  present  in  the  hydrates.  They  concluded  that 
there  is  some  evidence  of  0-H****X"  and  0-H****0  bonding. 

The  problem  of  whether  or  not  all  the  hydrogen  atoms  are 
Involved  in  hydrogen  bonding  was  called  to  our  attention  by 


Hornig's  private  communication,  during  the  summer  of  1961, 
to  G.  B.  Carpenter, 

By  examination  of  the  interatomic  distances  in 
the  sodium  bromide  dihydrate  crystal,  it  should  be  possible 
to  confirm  and  extend  the  findings  of  Schiffer  and  Hornig. 

It  was  also  desired  to  investigate  the  possibility  that, 
from  a  well-refined  crystal  structure  model  of  this  hydrate, 
the  hydrogen  atoms  may  be  evident  in  a  three  dimensional 
difference  electron  density  synthesis. 

During  the  final  stages  of  this  investigation,  it 
was  learned  that  Van  Meerssche,  Culot  and  Piret  (1962)  had 
recently  completed  their  independent  determination  of  the 
crystal  structure  of  sodium  bromide  dihydrate  by  X-ray  meth¬ 
ods.  The  reported  details  of  the  structure  are  essentially 
the  same  as  those  concluded  in  this  work.  More  detailed 
comparisons  will  appear  later. 


3. 


II »  Experimental  Procedures 

A.  Preparation  of  Sodium  Bromide  Dihydrate 

The  method  of  preparation  of  sodiuii  bromide  di¬ 
hydrate  (Gmelin,  1928,  p.  420-421)  is  simply  recrystal¬ 
lization  from  an  aqueous  solution  of  sodium  bromide, 

A  water  solution  of  sodivim  bromide,  saturated  at 
42“C,  was  prepared  in  a  beaker  which  was  kept  in  a  ther- 
mostatted  electric  oven.  By  controlled,  very  slow  cooling 
of  the  oven  down  to  room  temperature  (about  25“C),  large, 
irregular,  flat  plates  of  translucent  material  arranged  in 
step-wise  layers  crystallized  out  of  the  solution.  These 
plates  could  be  easily  broken  apart  by  mechanical  means  to 
yield  what  appeared  to  be  single  crystals.  The  percentage 
of  water  in  the  substance  was  determined  roughly  by  weighing 
a  sample  before  and  after  heating  in  an  oven.  It  corres¬ 
ponded  to  the  molecular  formula  NaBr*2H20. 

B.  Growth  and  Mounting  of  a  Single  Crystal. 

There  were  three  considerations  involved  in  the 
problem  of  growing  and  using  a  single  crystal  of  sodium 
bromide  dihydrate  in  an  X-ray  study.  Firstly,  if  a  crystal 
of  this  substance  were  left  exposed  to  the  air,  it  would 
either  deliquesce  or  effloresce  strongly,  depending  upon  the 
prevailing  atmospheric  conditions.  Secondly,  the  "optimum" 
cross-sectional  diameters  (assuming  a  cylindrical  crystal) 
are  0.142  mm  for  CuKoc  radiation  and  0.197  mm  for  MoKoc 
radiation.  (These  values  were  calculated  from  the  equation 
t  -  2/|a  (Buerger,  1942),  where  Is  the  linear  absorption 
ooef f Iclent  in  cm”^  calculated  from  ja.  *  d  ^  , 
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In  the  latter  equation,  d  is  the  density  of  the  crystal  in 
-3 

gm-cm  ,  Pj is  the  weight-fraction  of  the  ith  element  in  the 

compound,  and  )t  is  the  mass  absorption  coefficient  in 

-1  2 

gm  cm  ,  characteristic  of  the  ith  element  and  the  type  of  X- 
radiation  used.)  Growing  and  handling  of  a  very  small 
"optimum-size”  crystal  would  be  somewhat  awkward.  Finally, 
if  an  irregularly-shaped  crystal  were  used,  the  correction 
for  its  absorption  of  X-rays  would  be  complicated. 

These  problems  could  be  neatly  solved  by  growing 
single  crystals,  again  from  saturated  solutions  as  outlined 
above,  in  5-8  cm  lengths  of  thin-walled  Pyrex  capillary 
tubing  of  inside  diameter  about  0.1-0. 2  mm.  After  warm, 

saturated  solution  was  drawn  up  into  the  capillaries,  both 
ends  were  sealed  off  with  Picein  wax.  A  large  number  of 
these  capillaries  containing  crystalline  material  were  ex¬ 
amined,  both  visually  and  under  a  polarizing  microscope,  in 
order  to  detect  possibly  serious  defects  of  shape,  size,  or 
possible  twinning  in  the  crystals.  Finally,  a  practically 
cylindrical  crystal  measuring  about  0.19  mm.  in  diameter  and 
3  mm  in  length,  contained  in  a  5  cm  length  of  capillary, 
was  deemed  acceptable  for  the  subsequent  X-ray  photography. 
The  capillary  was  carefully  snapped  off  and  re-sealed  with 
Picein  wax  so  as  to  leave  a  3  cm  length  of  capillary,  with 
the  desired  crystal  about  1  cm  from  one  end.  The  other  end 
was  mounted  by  means  of  Picein  wax  on  the  end  of  an  8  mm . 
length  of  solid  glass  rod  3  mn  in  diameter.  This  rod  was 
then  inserted  in  the  opening  provided  on  a  goniometer  head. 
The  crystal  could  then  be  placed  in  the  path  of  the  X-ray 
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beam  by  screwing  the  whole  goniometer  head  Inte  place  on  the 
chosen  X-ray  diffraction  camera, 

C,  Recording  of  X-Ray  Diffraction  Spectra 
1,  Apparatus  and  Equipment  Involved 
a.  X-Ray  Sources 

The  X-rays  were  produced  from  sealed,  water- 
cooled  X-ray  tubes  set  in  a  Norelco  X-ray  Diffraction  Unit 
of  Type  12031,  manufactured  by  the  North  American  Philips 
Co.,  New  York,  N.Y.  Tubes  with  targets  of  molybednum  and 
of  copper  were  used.  Unfiltered  molybdenum  radiation  was 
employed  for  all  orientation  photographs,  but  practically 
monochromatic  radiation  was  desirable  for  other  purposes. 

MoKo<  radiation  was  partly  isolated  by  intercepting  the  X-ray 
beam  with  a  zirconium  filter,  CuKo<  by  interception  with  a 
nickel  filter.  The  molybdenum  radiation  was  produced  at  a 
filament  current  of  15-20  ma  and  voltage  of  55  kV,  the  copper 
radiation  at  20  ma  and  35  kV.  Molybdenum  radiation  was  used 
for  most  intensity  photographs  to  minimize  absorption  of  the 
X-rays;  copper  radiation  was  used  only  to  obtain  certain 
small-angle  reflections. 

b.  Cameras 

Use  was  made  of  the  Weissenberg  camera  and  of 
the  Buerger  precession  camera,  both  of  which  were  manufactured 
by  the  Charles  Supper  Co.,  Newton  Center,  Mass.  The 
Weissenberg  camera  accommodates  an  approximately  12i^  by 
14i  cm  rectangle  of  X-ray  film,  wJ  ^reas  12i  cm  squares  of 
X-ray  film  are  usually  fitted  into  the  precession  camera. 

In  practice,  the  film  which  is  loaded  into  the 
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V«issenberg  camera  is  translated  back  and  forth  as  the 
crystal,  bathed  in  the  X-ray  beam,  is  rotated.  (In  this 
case,  the  rotation  axis  was  the  long  axis  of  the  capillary 
containing  the  single  crystal  of  sodium  bromide  dihydrate.) 
By  use  of  a  cylindrical  metal  screen  called  a  layer-line 
screen,  it  is  possible  to  limit  the  film-recorded  reflec¬ 
tions  to  those  in  a  given  plane  of  the  reciprocal  lattice. 
The  resulting  developed  photograph  shows  "festoons”  of 
spots  which  represent  a  distorted  picture  of  a  single 
plane  of  the  reciprocal  lattice.  The  distortion,  however, 
does  not  prevent  one  from  assigning  hkl  indices  to  each 
spot  with  relative  ease.  A  comprehensive,  detailed  ex¬ 
planation  of  the  principles  and  geometry  of  the  Weissenberg 
camera,  the  recording  of  reflections  by  the  equi-inclination 
method,  and  the  simple  technique  for  rapid  indexing  is 
presented  by  Buerger  (1942,  pp.  214-295). 

The  geometry  of  motion  involved  in  the  preces¬ 
sion  method  of  photography  is  somewhat  complicated;  the 
normal  to  the  reciprocal  lattice  plane  being  photographed 
sweeps  out  a  circular  cone  about  the  X-ray  beam.  The  pre¬ 
cession  method  has  the  advantage  of  providing  undistorted 
pictures  of  the  reciprocal  lattice,  thus  facilitating  the 
rapid  and  accurate  determination  of  unit  cell  constants, 
as  well  as  making  the  indexing  quite  obvious.  A  further 
advantage  is  seen  in  that  it  is  possible  to  photograph 
several  groups  of  planes  while  the  crystal  is  mounted  in  a 
single  orientation.  There  is,  however,  a  blind  spot 
present  at  the  center  of  the  diffraction  pattern  of  each 
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level,  and  this  blind  spot  Increases  rapidly  in  area  as  the 
distance  from  the  reciprocal  lattice  origin  increases.  Full 
details  of  the  precession  method  are  presented  in  the  Supper 
Co.  Manual  and  in  an  ASXRED  monograph  (Buerger,  1944). 

Both  Weissenberg  and  precession  intensity  photographs  were 
taken  in  the  course  of  this  study. 

c.  Film.  Developing  Materials  and  Procedure 
No-Screen  Medical  X-Ray  Film,  manufactured  by 
Eastman  Kodak  Co.,  Rochester,  N.Y.,  was  used  for  all  X-ray 
photographs.  Black  cellulose  acetate  plastic  was  used  to 
shield  the  X-ray  film  from  light  and  to  hold  the  film  in 
place  in  the  camera. 

The  developing  procedure  for  all  photographs 
was  as  follows:  The  exposed  film  was  immersed  for  5i  minutes 
in  Kodak  X-Ray  Developer  treated  with  Kodak  AntiFog  No.  2 
(0.3  g/l  of  developer),  held  in  a  dilute  acetic  acid  "stop” 
bath  for  about  3  seconds,  immersed  in  Kodak  X-Ray  Fixer  for 
twenty  minutes,  and  washed  finally  in  cold  running  water  for 
30  minutes. 

2.  X-Ray  Diffraction  Photographs 
a.  Weissenberg  Photographs 
The  crystal  was  oriented  on  the  Weissenberg 
apparatus  by  making  use  of  20"  oscillation  photographs  taken 
90®  apart.  By  adaptation  of  Bunn's  method  (1946,  p. 173-174 ) , 
it  was  possible  to  estimate  the  correction  to  be  applied  to 
the  crystal  setting  in  order  that  the  rotation  axis  and  the 
£  (crystallographic)  axis  should  coincide.  (The  c  axis  was 
so  designated  to  facilitate  comparison  between  the  cell 


parameters  obtained  in  this  study  and  those  obtained  by 
Wooster  (1932)).  The  orientation  of  the  crystal  was 
checked  by  examination  of  a  6-hour  rotation  photograph. 

Two  series  of  24-hour  Weissenberg  intensity  photographs  of 
the  a*-b*  plane  were  made  with  an  oscillation  ranee  of  about 
200“.  The  first  series  was  made  using  MoKoc  radiation; 
levels  0,1, 2, 3, 4,  and  5  were  included  in  this  series.  The 
second  series  was  made  using  CuKoj  radiation;  levels 
0,1,2,  and  3  were  included  in  this  series.  These  two  sets 
of  photographs  complemented  each  other  inasmuch  as  the  MoKoc 
set  provided  good,  large  angle  intensity  data,  while  the 
CuKoc  set  provided  small  angle  reflections  not  available 
on  the  MoKo,  set.  Both  series  employed  the  multiple  film 
technique  as  described,  for  example,  by  Richards  (1959). 

For  MoKck  radiation,  the  pack  of  three  films  was  separated 
by  two  thin  sheets  of  brass  foil  before  being  loaded  into 
the  camera.  But  for  both  MoKo?  and  CuKoc  radiation,  the 
ratios  of  the  intensities  of  a  given  diffraction  spot  re¬ 
corded  simultaneously  on  the  three  films  was  about  16:4:1. 

In  this  manner,  most  of  the  weakest  and  strongest  intensi¬ 
ties  were  conveniently  measurable.  All  of  the  intensities 
belonging  to  one  given  reciprocal  lattice  level  were  sub¬ 
sequently  placed  upon  the  same  numerical  scale  by  the  appli¬ 
cation  of  film  scale  factors. 

b.  Precession  Photographs 
The  crystal  was  transferred  to  the  precession 
apparatus  for  orientation.  Since  the  crystal  had  already 
been  oriented  for  Weissenberg  photography,  only  the  dial 
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axis  of  the  precession  instrument  needed  adjustment  for  the 
a*-c*  plane  to  be  in  photographic  position.  The  amount  of 
this  adjustment  could  be  found  by  examination  of  a  zero- 
level  Weissenberg  exposure  for  the  position  of  reflection  of 
the  spots  aligned  along  the  b*  axis.  Final  adjustments  were 
carried  out  with  the  aid  of  the  "Orientation  Error  Chart  for 
Precession  Method"  (Supper  Go.  Ma  lual).  Measurements  ob¬ 
tained  from  a  cone-axis  photograph,  coupled  with  the  ap¬ 
plication  of  a  standard  chart  (Evans,  Tilden  and  Adams, 

1949),  provided  camera  settings  which  led  to  the  precession 
photography  of  five  reciprocal  lattice  levels — levels 
0,1, 2, 3  and  4.  Exposures  of  24  hours,  6  hours  and  li  hours 
were  made  of  each  level.  It  should  be  noted  that  about  mid¬ 
way  throu^  the  precession  intensity  series,  a  whitish  cast 
had  developed  upon  the  surface  of  the  crystal,  indicating 
that  some  deterioration  of  the  crystal  had  taken  place 
during  its  exposure  to  X-rays.  That  the  deterioration  did 
not  significantly  affect  the  diffraction  properties  of  the 
crystal  was  shown  by  the  comparison  of  a  first-level  photo¬ 
graph  taken  before  the  whiteness  was  noticed  with  one  taken 
well  after  it  became  visible. 

D.  Determination  of  Unit  Cell  Dimensions 
1.  Film  Measurements 

Preliminary  examination  of  the  Weissenberg  and 
precession  records  showed  that  Wooster’s  (1932)  monoclinic 
unit  cell  is  correct.  His  choice  of  b  as  the  unique  axis 
was  retained. 

Measurements  were  performed  on  a  zero-level 
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precession  film  in  order  to  determine  a*  £*  and  cosine  ^  *. 
Another  zero-level  precession  film,  obtained  after  rotating 
the  precession  spindle  axis  90®,  was  used  to  determine  b*’'. 
Seventy-two  spot  positions  and  twelve  spot  positions,  respec¬ 
tively,  were  measured  to  0.01  mm  v;ith  the  aid  of  a  film 
viewer  manufactured  by  the  Charles  Supper  Co.,  Newton 
Center,  Mass. 

2.  Calculating  the  Cell  Dimensions 

The  reciprocal  cell  parameters  were  calculated  by 
application  of  the  method  of  least  squares,  vdiich,  in  this 
case,  minimized  the  quantity  In  this 

expression,  Sg  is  the  observed  squared  reciprocal  lattice 
vector  calculated  for  reflection  hkl  by  the  quantity 
4  sin2  0/;^»  ,  where  Q  is  the  Bragg  angle  for  reflection 

hkl,  and  ^  is  the  wavelength  of  X-radiation  used  in  record¬ 
ing  the  reflection;  wjjki  is  the  v/eight  assigned  to  S^.  The 
sum  is  taken  over  all  the  measured  reflections.  Weights  were 
subjectively  assigned  on  the  basis  of  the  magnitude  of  S  and 
ease  of  measurement  of  the  spot  position.  It  was  a  simple 
matter  to  convert  the  reciprocal  cell  constants  into  those 
of  the  direct  lattice. 

The  standard  deviations  associated  with  ^ ,  £*  and 
^  *  were  calculated  by  means  of  the  equation  for  the  vari¬ 
ance,  <ri*  ,  which  is 

(Whittaker  and  Robinson,  1940),  where  m  is  the  number  of  ob¬ 
servational  equations,  s  is  the  nixmber  of  parameters  involved, 
is  the  ith  cofactor  of  the  determinant  D  of  the  coef¬ 
ficients  of  the  normal  equations,  and  (^So.  —  is  ^he 
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residual  of  the  jth  observational  equation.  The  standard 
deviation  of  b*  was  obtained  in  like  manner.  Again, 
the  standard  deviations  in  the  direct  cell  constants  may 
be  easily  obtained  from  these.  It  must  be  realized  that 
the  standard  deviations  in  the  cell  constants  are  measures 
only  of  the  internal  consistency  of  the  data,  and  do  not 
take  systematic  errors  into  account.  (The  sources  of 
such  errors  may  be,  for  example,  film  shrinkage,  or 
failure  to  calibrate  the  X-ray  cameras.  These  factors 
were  not  corrected  for  in  this  work.) 

3.  Discussion  of  the  Cell  Dimensions 

The  only  X-ray  work  on  sodium  bromide  dihydrate 
published  prior  to  the  beginning  of  this  study  was  a  deter¬ 
mination  of  the  cell  constants  and  space  group  by  Wooster 
(1932).  He  employed  Weissenberg  and  oscillation  photo¬ 
graphs  for  this  purpose. 

Very  recently.  Van  Meeresche,  Culot,  and  Piret 
(1962)  have  also  determined  cell  constants  (in  the  course 
of  their  complete  structure  determination)  by  making  use 
of  Weissenberg  and  rotation  photographs.  For  comparison. 
Table  I  below  lists  the  cell  parameters,  number  of 
molecules  per  unit  cell,  and  the  calculated  densities 
resulting  from  the  three  independent  investigations.  Note 
that  the  limits  given  in  this  work  are  standard  devia¬ 
tions.  The  actual  errors  are  probably  of  the  same  order 
of  magnitude  as  those  of  Van  Meersche  and  co-workers. 
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Upon  examination  of  Table  I,  it  is  clear  that 
Wooster's  set  of  cell  constants  do  not  approach  the  value 
"4"  for  the  calculated  number  of  molecules  per  unit  cell 
quite  as  well  as  do  the  other  two  sets.  Moreover,  the 
agreement  is  quite  good  between  the  individual  axial  lengths 
and  monoclinic  angle  given  here  and  those  given  by  Van 
Meerssche  and  co-workers.  For  these  reasons,  it  would  seem 
that  either  of  these  two  sets  of  cell  parameters  are 
preferable  to  Wooster's. 

TABLE  I 

Comparison  of  Unit  Cell  Quantities 
Quantity  This  Thesis  Work  Van  Meerssche  Wooster 

(!C96^  and  co-workers  ( 1932 ) 

(1962) 

a  6.575±0.002A  6.57±0.03A  6.59J?(  ±0.5?S) 

b  10.456±0.005X  10.38t0.05X  10.20X{±0.5?9 

c  6.776t0.002X  6.78±0.03X  6.5lA(±0.5?C) 

(3  113‘»23'±0.002“  113»30'±0.4“  112'»5'  (tO.SjS) 

Calc'd  No. 

Molecules 

per  Unit  Cell  4.04  4.00  3.83 

Calc'd  23-25®C  unspec'd  unspec'd 

Density  2.16  room  temp.  2.18  room  temp.  2.16  ■'oom  temp. 

Experimental  Density  is  D|0  =  2.176t0.003  (Gmelin,  1928) 
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E.  Determination  of  the  Space  Group 

It  was  already  indicated  that  the  cell  parameters 
of  sodium  bromide  dihydrate  are  characteristic  of  the  mono¬ 
clinic  crystal  class.  In  addition,  it  was  noted  upon  ex¬ 
amination  of  the  X-ray  photographs  that^  for  the  class  of  hOl 
reflections,  only  those  for  vAiich  l»2n  (n  is  an  integer) 
were  present,  and  that  for  the  class  of  OkO  reflections,  only 
those  for  which  k  =  2n  were  present.  It  was  then  possible  to 
determine  the  space  group  corresponding  to  these  systematic 
absences  as  P2i/c  uniquely  (Internationale  Tabellen,  1935, 
p.  382).  This  finding  is  in  agreement  with  those  of 
Wooster  and  of  Van  Meerssche. 

F.  Obtaining  Structure  Factor  Magnitudes  from  Corrected 
Relative  Intensity  Data 

1.  Microphotometer  Method  for  Measuring  Intensities 
on  Precession  Films 

All  intensity  measurements  of  diffraction  spots  on 
precession  films  were  made  with  the  aid  of  a  Recording  Micro¬ 
photometer  manufactured  by  the  F.  C.  Henson  Co.,  Pasadena, 
California,  'fhe  original  recording  and  detection  system  was 
modified  to  consist  essentially  of  a  light  source  and  a 
photomultiplier  tube  electrically  connected  to  a  micro- 
ammeter.  The  procedure  of  estimating  an  intensity  was  as 
follows:  The  spot,  the  intensity  of  which  v^as  desired,  was 

positioned  so  that  the  area  of  the  spot's  maximum  blackness 
intercepted  a  light  beam  which  impinged  upon  a  photoelectric 
cell.  The  amount  of  light  viiich  passed  through  the  spot 
was  indicated  by  the  deflection  of  a  sensitive  microammeter. 
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The  background  (that  is,  the  intensity  of  the  area  im¬ 
mediately  surrounding  a  spot)  was  also  measured  in  the  same 
way.  The  relative  intensity  of  the  given  spot  is  then 
Irel  ■  log  (average  galvanometer  reading  corresponding  to 
background)  -  log (galvanometer  reading  corresponding  to 
spot).  In  this  study,  the  differences  were  multiplied  by 
100  for  numerical  convenience.  It  was  necessary  to  measure 
only  half  of  the  reflections  recorded  on  the  intensity 
photograph;  one  half  of  a  given  film  was  related  to  the  other 
half  by  a  center  of  symmetry. 

2.  Visual  Method  for  Measuring  Intensities  on 
Weissenberg  Films 

An  adaptation  of  the  visual  method  described,  for 
example,  by  Lee  (1958)  was  used  in  this  study,  since  the  ar- 
rargement  and  shape  of  the  Weissenberg  spots  did  not 
facilitate  the  measurement  of  their  intensities  by  the 
microphoto meter  method.  In  addition,  the  opportunity  of 
comparing  the  intensities  resulting  from  the  two  different 
methods  of  measurement  is  presented. 

The  visual  method  consists  in  comparing  the 
intensity  of  a  given  spot  to  the  intensities  of  one  or  two 
spots  belonging  to  a  series  of  spots  of  assigned,  graded, 
relative  intensities.  By  interpolation,  an  estimated 
value  of  relative  intensity  is  assigned  to  the  given  spot. 

In  order  to  obtain  the  desired  series  of  spots  for  pur¬ 
poses  of  comparison,  it  was  necessary  to  prepare  an 
••intensity  strip*’.  The  intensity  strip  was  made  as 


follows: 


A  region  of  reflection  within  the  crystal  cor¬ 
responding  to  a  suitable  diffraction  spot  on  the  film 
(suitable  in  blackness,  absence  of  strong  background,  and 
shape)  was  isolated  within  a  4-8“  oscillation  on  the  Weis- 
senberg  apparatus.  The  reflection  region  was  exposed  to 
X-rays  for  a  pre -determined  length  of  time  (or  number  of 
oscillations)  causing  the  corresponding  spot  on  the 
processed  film  to  be  blackened  accordingly.  The  film  was 
given  a  3  mm  translation,  end  the  reflection  region  again 
exposed  to  X-rays  for  another  known  length  of  time.  The 
process  was  repeated  until  a  series  of  20  spots,  arranged 
in  a  strait t  line  and  graded  in  blackness,  were  produced. 
The  time  intervals  of  X-ray  exposure  followed  a  geometric 
series  of  form  (base)^  ■  exposure  time,  where,  in  this 
study,  the  base  is  1.25,  x  assumes  all  integral  values  from 
0-19,  and  exposure  time  is  in  minutes.  In  principle,  any 
numerical  base  with  any  regular  sequence  of  exponents  could 
be  utilized,  and  the  result  could  be  expressed  in  any  time 
units,  depending  upon  the  nature  of  the  crystal  used.  A 
geometric  series  of  exposure  times  was  employed  because  the 
neuro-physiology  of  the  eye  is  such  that  it  responds  to  the 
logarithm  of  the  exposure  time  rather  than  to  the  exposure 
time  itself  when  estimating  intensities.  Thus,  if  r-  geo¬ 
metric  series  of  time  exposures  is  chosen,  the  eye  sees  the 
same  interval  of  blackness  between  any  two  successive  spots 
in  the  intensity  strip. 

Usually,  vhen  the  relative  intensity  of  a  given 
spot  was  estimated,  it  was  between  two  spots  (in  blackness) 


16 


on  the  intensity  strip.  To  obtain  the  desired  relative 
intensity,  use  was  made  of  the  relation  log  Ix“r  los  Id 
(1-r)  log  Ii,  where  is  the  desired  relative  intensity  of 
the  given  spot,  is  the  relative  intensity  of  the  darker 
spot,  Ii  is  the  relative  intensity  of  the  lighter  spot,  and 
r  is  the  fraction  of  interpolation.  It  was  found  convenient 
to  construct  a  table  of  values  of  Ix  for  the  r  values  of 
1/5,  1/4,  1/3,  2/5,  1/2,  3/5,  2/3,  3/4  and  4/5.  Thus,  all 
that  needed  to  be  done,  in  pracuitfer,  was  to  write  down  the 
exponent  corresponding  to  the  lighter  reflection  together 
with  the  fraction  of  interpolation,  and  Ix  could  then  be 
easily  read  off  from  the  table.  The  intensity  of  every 
spot  that  could  be  detected  by  eye  on  the  Weissenberg  films 
was  visually  estimated. 

All  Weissenberg  and  precession  intensity  measure¬ 
ments  were  performed  in  no  particular  order,  and  under  con¬ 
ditions  which  were  as  uniform  as  possible.  Film  factors  were 
computed  after  all  intensity  measurements  were  completed  in 
order  to  eliminate  bias. 

3.  Comparison  of  the  Micro photometer  and  Visual 

Methods 

The  micro photometer  method  yields  measurements 
of  slightly  greater  precision  than  are  obtained  by  the 
visual  method.  A  micro photometer  measurement  may  enjoy 
5-10^  reproducibility  whereas  a  visual  measurement  ranges 
from  5-15^  in  reproducibility.  Furthermore,  the  micro¬ 
photometer  method  depends  less  upon  the  physical  condition 
of  the  experimenter  and  less  upon  local  lighting  conditions. 
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On  the  other  hand,  it  must  be  noted  that  the  micro¬ 
photometer  method  yields  peak  intensities  rather  than  the  more 
desirable  integrated  intensities  given  by  the  visual  method. 

In  addition,  it  cannot  be  denied  that  the  visual  method  is 
faster,  enabling  one  to  perform  2  or  3  measurements  in  the 
same  time  required  for  one  measurement  by  the  micro photometer 
method.  This  fact  is  of  some  importance  if  one  is  interested 
in  obtaining  the  intensities  of  a  thousand  or  more  diffraction 
spots.  Finally,  serious  errors  are  more  likely  to  be  avoided 
if  one  is  able  to  measure  an  hkl  reflection  several  times  by 
the  visual  method  rather  than  once  or  twice  by  the  slightly 
more  precise  microphotometer  method. 

Taking  all  factors  into  consideration,  it  is  the 
author’s  opinion  that,  in  spite  of  the  slight  superiority  in 
precision  of  the  micro photometer  method,  and  that  provided 
good  intensity  strips  can  be  readily  made,  the  visual  method, 
which  produces  rapid  estimates  of  integrated  intensities,  is 
generally  the  more  desirable  method. 

4.  Extraction  of  Structure  Factor  Magnitudes  from 
Relative  Intensities 

Squared  structure  f.actor  magnitudes  were  ob¬ 
tained  by  multiplying  the  relative  intensities,  the  Ihkl's, 
by  the  reciprocals  of  the  Lorentz  and  the  polarization 
factors.  Mathematically,  iFhkll  ^  ■  (^p)  jj^l^lhkl  (Buerger, 

1960,  pp.  25-48),  where  iFy^kll  is  the  magnitude  of  the  struc¬ 
ture  factor  corresponding  to  reflection  hkl,  and  (Lp)hjcl 
represents  the  reciprocal  product  of  the  Lorentz  factor  (in¬ 
cluding  any  geometrical  factors  associated  with  the 
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particular  method  of  X-ray  photography  employed)  and  the 
polarization  factor. 

A  set  of  Lorent z-polarization  charts  reproduced 
upon  transparent  film  were  available  for  use  in  applying  the 
Lorentz-polari zation  correction  tr  the  relative  intensity 
data  obtained  frcSHrprecession  films.  Save  for  adjustments 
of  scale,  these  charts  are  identical  to  those  made  by  Waser 
(1951),  and  by  Grenville-Wells  and  Abrahams  (1952).  Waser 
constructed  a  Lorentz-polarization  chart  for  the  zero  level 
of  the  reciprocal  lattice  and  precession  angle  =  .’30®; 
the  seven  Grenville-Wells  charts  are  constructed  at  inter¬ 
vals  of  0.05  reciprocal  lattice  units  above  the  zero 
level,  again  for  ~  30°.  Linear  interpolation  was  applied 
to  obtain  Lorentz-polarization  factors  for  precession  films 
corresponding  to  intermediate  reciprocal  lattice  levels. 
Extraction  of  the  square  roots  of  the  squared  stznicture 
factor  magnitudes  resulted  in  a  set  of  "precession”  struc¬ 
ture  factor  magnitudes. 

The  Lorentz-polarization  corrections  for  rela¬ 
tive  intensities  from  Weissenberg  films  were  applied  by  com¬ 
puter  means  at  the  Brown  University  Computing  Laboratory.  A 
program  developed  for  this  purpose  was  utilized  (Carpenter 
and  Hall,  1962,  pp.  8-12).  A  set  of  "Weissenberg"  structure 
factor  magnitudes  along  with  their  corresponding  squared 
values  were  obtained  as  computer  output. 

5.  Placing  All  Structure  Factor  Magnitudes  on  the 

Same  Scale 


In  the  preceding,  it  was  explained  how  relative 
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structure  factor  magnitudes  were  obtained  for  reflections 
appearing  on  films  corresponding  to  seven  reciprocal  lattice 
levels  photographed  by  the  Weissenberg  method,  and  for  re¬ 
flections  appearing  on  films  corresponding  to  five  reciprocal 
lattice  levels  photographed  by  the  precession  method. 

The  next  problem  was  that  the  structure  factor 
magnitudes  corresponding  to  one  film  were  on  a  numerical 
scale  different  from  that  for  any  other  of  the  twelve 
films,  It  was  necessary  to  place  all  the  structure  factor 
magnitudes  on  the  same  numerical  scale  by  finding  and  ap¬ 
plying  twelve  scale  factors,  k]^,  k2 ,  k3, . ^12  > 

structure  factor  magnitudes  corresponding  to  the  twelve 
films.  Evaluation  of  these  scale  factors  depended  on  the 
fact  that  a  considerable  number  of  reflections  were  re¬ 
corded  at  two  film  positions,  once  on  a  precession  film  and 
again  on  a  Weissenberg  film. 

In  order  to  find  the  set  of  k’s  which  would  lead 
to  the  best  adjustment  of  all  the  structure  factor  magni¬ 
tudes,  it  was  decided  to  apply  a  modification  of  the  method 
of  Rollett  and  Sparks  (1960).  The  finding  of  the  best  set 
of  k’s  became  essentially  a  least-squares  problem.  The 
function  v;hich  was  minimized  with  respect  to  the  k’s  is 

.  r  ^  •  Sub¬ 

scripts  i  and  j  refer  to  two  different  films  of  a  total  of 
N  different  films  to  be  correlated;  Sij  designates  reflec¬ 
tions  common  to  films  i  and  j;  Fi(Sij)  is  the  magnitude  of 
the  structure  factor  for  reflection  as  measured  on  film 
i;  ki  is  the  scale  factor  to  be  applied  to  the  structure 
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factor  magnitudes  of  film  i;  w(Sj^j)  is  the  weighting  factor 
to  be  applied  to  reflection  It  wa.s  necessary  to  intro¬ 

duce  an  undetermined  multiplier  to  find  a  set  of  k*s  other 
than  ki  ■  0  for  all  i.  The  problem  then  became  one  of  solving 
a  12X12  maurix  ^or  the  minimum  eigenvalue  and  its  associated 
eigenvector  (or  set  of  scale  factors).  Fortimately ,  the 
Brown  Computing  Laboratory  facilities  together  with  an 
eigenvalue-eigenvector  program  written  by  T.  Hughes  were 
available.  Thus,  after  the  proper  matrix  coefficients  were 
set  up  with  aid  of  a  desk  calculator,  the  eigenvalue  and 
eigenvector  were  readily  found. 

After  application  of  the  scale  factors  to  the 
films,  it  became  evident  that  certain  corrections  should 
have  been  applied  to  the  "raw”  (that  is,  uncorrected  for  the 
Lorentz-polarization  factors)  intensity  data  of  the  pre¬ 
cession  films.  Since  only  peak  precession  intensities  had 
been  measured  by  the  microphotometer  method,  it  was  neces¬ 
sary  to  correct  for  the  area  of  the  spot.  (The  need  for 
this  correction,  however,  does  not  arise  if  one  measures  in¬ 
tensities  by  the  visual  method,  owing  to  the  eye’s  natural 
allowance  fcr  area.)  The  correction  was  carried  out  by  the 
following  semi-empirical  means. 

Assume  a  roughly  conical  distribution  in  the 
diffraction  spot  intensity.  The  microphotometer  method 
yields  a  measure  of  the  cone  height,  h,  only.  Now  consider 
the  integrated  intensity  of  a  diffraction  spot  to  be 
rou^ly  analogous  to  the  volume  of  a  cone.  In  other  words, 
assume  that  the  raw,  integrated  intensity,  I,  is  given  by  the 
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product  of  the  height,  h,  and  base  area,  A,  of  the  cone. 

But  since  A  may  be  expressed  empirically  in  terms  of  h  by 
the  relation  h  ■  A^,  =  hA  =  hence  n  may  be 

evaluated  by  a  log  h  vs.  log  A  plot.  About  12-15  pre¬ 
cession  spots,  varying  in  blackness,  were  chosen  for  the 
plot.  The  spot  areas  were  measured  with  a  travelling 
microscope.  The  log-log  plot  showed  considerable  scatter, 
but,  nevertheless,  it  was  possible  to  draw  wha^t  appeared  to 
be  a  reasonable  straight  line  through  the  distribution  of 
points.  This  line  led  to  the  expression  for  the  corrected 
raw  intensity  which  was  used  to  correct  all  the 

previously-measured  precession  intensities.  A  new  set  of 
structure  factor  magnitudes  was  derived  from  the  corrected 
precession  intensities.  After  the  eigenvalue-eigenvector 
program  was  applied,  using  the  corrected  precession  struc¬ 
ture  factor  magnitudes,  a  set  of  k’s  resulted  which  led  to 
reasonably  well-correlated  structure  factor  magnitudes.  In 
this  way,  structure  factor  magnitudes,  all  on  one  numerical 
scale,  were  obtained  for  634  independent  reflections.  At 
this  point,  it  was  possible  to  begin  the  detailed  structure 
investigation  of  sodium  bromide  dihydrate.  . 
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III .  Determining  the  Detailed  Structure  of  Sodium 

Bromide  Hydrate 

A.  The  Phase  Problem  of  X-Ray  Crystallography 

The  structure  factor,  Fhkl>  in  general,  a 
complex  quantity  which  may  be  expressed  as  a  product  of 
(Fhkllj  ■'"he  structure  factor  magnitude  for  reflection  hkl, 
and  e^^hkl,  the  phase  factor.  If  the  complex  structure 
factor,  F^kli  were  experimentally  available,  solving  a 
crystal  structure  would  be  a  routine  matter.  Knowledge 
of  the  Fhkl's  would  enable  one  to  evaluate  the  electron 
density  P(x,y,z),  at  any  point  x,y,z  in  the  unit  cell  from 
the  relation  ® '/y 

Cochran,  1955,  p.  12),  the  summation  extending  over  all  in¬ 
dices  hkl.  An  electron  density  map  could  tiien  be  constructed, 
which  would  reveal  the  atomic  positions  in  the  unit  cell. 

But  thus  far,  it  has  been  possible  to  obtain  only  structure 
factor  magnitudes  by  experimental  means.  It  would  seem  that, 
owing  to  the  lack  of  phase  information,  the  solution  of  a 
crystal  structure  would  be  at  least  very  difficult.  This 
constitutes  the  phase  problem  of  X-ray  crystallography  ex¬ 
pressed,  for  example,  by  Buerger  (1960,  p.  551).  One 
attack  on  the  phase  problem  is  by  means  of  the  Patterson 
(or  \F)2)  synthesis. 

B.  The  Patterson  Svnthesi.s 

1.  The  Physical  Significance  of  the  Patterson  Function 
An  insight  into  the  physical  significance  of  the 
Patterson  function,  P(u,v,w) ,  is  provided  by  examination  of 
one  expression  of  P(u,v,w)  which  is 
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P ( x+u , y+v , z+w ) dxdydz . 

This  relationship  expresses  the  fact  that  the  Patterson 
function,  evaluated  for  the  chosen  vectorial  distance 
u,v,w,  is  given  by  the  product  of  i^(x,y,2),  the  electron 
density  at  x,y,z  and  P (x+u, y+v , z+w) ,  the  electron  density 
at  x+u, y+v, z+w.  Fractional  coordinates  x,y,z  may  vary 
through  all  values  from  0  to  1,  v*iile  u,v,w  are  held  fixed. 

The  highest  positive  numerical  value  of  the  Patterson  function 
will  be  at  those  particular  values  of  u,v,  and  w  which  are 
coordinate  differences  between  peaks  in  the  electron  density, 
or  which,  in  other  words,  correspond  to  certain  interatomic 
vectors.  Whereas  an  interatomic  vector  may  have  its  origin 
anywhere  in  the  unit  cell,  that  vector  must  begin  at  the 
origin  in  "Patterson  (or  vector)  space".  The  preceding  was 
condensed  from  Buerger's  book  (195'^,  chapter  2),  which  con¬ 
tains  a  detailed  discussion  of  the  theory  and  application  of 
the  Patterson  function. 

2.  Forms  of  the  Patterson  Function  Suitable  for 
Calculation 

The  general  expression  for  the  Patterson  function 
is  given  as  (Lipson  and  Cochran,  1953,  pp.  9-15) 

P(u,v,w)  =  l/vy  /  /  |Fhkl)^cos2Tf(hu+kv+lw) ,  where  u,v,w 
refer  to  fractions  of  the  lengths  of  a,b,c  axes,  respective¬ 
ly,  V  is  the  unit  cell  volume,  and  the  triple  summation  is 
taken  over  indices  h,k,  and  1.  By  taking  advantage  of  the 
space  group  symmetry  of  the  particular  substance  under  study, 
the  general  expression  is  usually  modified  to  facilitate 
practical  calculation.  The  modifications  involved  in  the 
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calculation  of  the  Patterson  projection  and  the  three- 
dimensional  Patterson  synthesis,  both  for  space  group 
P2j^/c ,  vdll  be  considered  next. 

Before  all  the  intensity  data  had  been  gathered,  a 
zero-level  Patterson  projection  was  constructed,  using  the 
available  squared  structure  factor  magnitudes  of  the  zero- 
level  precession  film,  in  order  to  obtain  preliminary  in¬ 
dications  of  the  structure  of  sodium  bromide  dihydrate. 

Since  a  Fourier  synthesis  computer  program  had  not  been 
completed  at  that  time,  the  Patterson  projection  was  cal¬ 
culated  with  Beevers-Lipson  strips  (Lipson  and  Cochran, 

1953,  pp.  39-93).  The  general  form  of  the  Patterson 
function  (as  given  above)  had  to  be  modified  in  the  following 
manner  to  allow  the  efficient  use  of  this  strip  method. 

Since  the  projection  was  along  the  b  axis  onto  the  zero  level, 
k  »  0,  and  the  general  form  could  be  written  as 
Pej/A  /  |F]^q2]  ^  cosBTT (hu+lw) .  (This  form  may  be 
derived  by  observing  that  an  electron  density  projected  in 
a  certain  direction  corresponds  mathematically  to  an  inte¬ 
gration  over  the  lattice  period  in  that  direction  (Buerger, 
1960,  p.  382)).  In  this  case,  ?[x,z)  P {x,y,z)bdy,  where 
b  refers  to  the  axis  along  which  the  projection  is  taken, 
P(x,z)  is  th#  projected  electron  density  at  point  x,z,  and 
P(x,y,z)  is  the  electron  density  at  x,y,z.  The  form  of 
Patterson  projection  oould  then  be  obtained  by  a  treatment 
similar  to  that  given  by  Lipson  and  Cochran  (1953,  pp. 

11-13).)  Application  of  the  trigonometric  identity 
cos  ( ^  )  -  cos  Of  cos^  -sin®C  sin^  resulted  in 


25 


the  form 

^  ,  |Fy^nil^^os(2'Trhu)cor.  (2Trlw)“Sin(2Tf  hu)  sin(2T(  Iw)]. 
It  was  necessary  to  measure  the  intensities  appearing  on 
only  one-half  of  the  film  (see  section  IIF  1),  hence  the 
summation  could  be  split  up  to  appear  as  1/4F  = 

O* 

[_cos(2Trhu)cos(2'n'lw)-sin{2'TT  hu)sin(2Tf  Iw)]  + 
1/A  ^os(2Tt  hu)  cos  ( 2Tf  Iw)  +sin{2Tt  hu)  sin(  2Tf  lw)J  . 

By  rearrangement  and  collection  of  terms,  the  Patterson 
function  for^the  zero-level  projection  could  be  written  as 
F  =  4/A  l^ROll  ^]cos(2'TThu)cos(2T(  Iw)  - 

[ifhOll  ^  -  IfhOll  sin(2T(hu)sin(2Tnw)^  .  Thus  it  was 
possible  to  evaluate  numerically  the  Patterson  function  at 
intervals  of  1/60  in  u  and  w  by  means  of  the  Beevers- 
Lipson  strips.  A  Patterson  map  was  then  constructed  by 
drawing  contour  lines  which  connected  points  of  equal 
"heights”. 

The  estimation  and  subsequent  processing  of  all 
of  the  spot  intensities  were  completed  at  about  the  same 
time  that  the  BXFS  computer  program  was  finished  (Carpenter 
and  Wilmot,  1962).  It  provides  a  rapid  means  of  evalu¬ 
ating  the  crystallographic  functions  usually  expressed  in 
terms  of  the  Fourier  series,  such  as  the  Patterson, 
electron  density,  difference  Patterson,  and  difference 
electron  dens5.ty  functions.  In  the  program,  the 
Fourier  synthesis  is  evaluated  in  the  form 

i^nBTrhxj^gJgzTTkyj^IJSarrl^,  where  G(hkl) 
is  a  linear  combination  of  either  structure  factors  or 
difference* between  structure  factors,  or  squared  structure 


26. 


factor  magnitudes,  or  differences  between  squared  struc¬ 
ture  factor  magnitudes.  In  order  to  obtain  the  Patterson 
function  in  this  form,  the  general  form  was  transformed  as 
will  now  be  explained. 

For  space  group  PE^/c,  where  k+l=2n  (n  being  an 

integer),  ^hkl"  ^hkl  ^^kl“  ^kl’  where 

k+l=2n+l,  again  F^kl=Fhkl  and  F^kl  “  F^kl  (International 
Tables,  1952,  p.  583).  Obviously,  it  was  sufficient  to 
deal  only  with  the  case  for  which  k+l=2n.  In  this  case 
the  summation  may  be  split  and  the  trigonometric  identity 
cos  (oc  +  ^  )  =^cos®f'cosP -sin*t^sinP  applied  to  yield 
P»l/2V^|^'^^'^j^Fhkltcos2Tt  {hu+lw)cos2Ttkv  - 
sin2TT(hu+lw)sin2Tf  kv]]  + 

F^l^os2TT  (hu+lw)cos2Tr(  -kv)-sin2TY{hu+lw)  sin2  Tt  (-ky)3  + 
FfJkl[cos  2Tt( -hu+lw)cos2Trkv  -  sin2TT(-hu+lw)sin2'Trkv'}  + 

-  sinETt  (hu-lw)sin2‘nk'0|  . 

After  combining  terms,  and  noting  that  the  case  k=2n+l 
yields  the  same  results,  the  expression  for  the  Patterson 
function  could  be  written  as 

P“4/V^^  ^  .^(Phyi  cos2Tf  (hu+lw)+F|klcos2Tf  (hu-lw)\cos2Tf  kv. 
o  * 

When  the  trigonometric  identities  cos(°f  +  ^)  =  cosOfcos^  - 
sin®<sin(?  and  cos(o<-^  )  =  cosC<cos^  +  sinorsin^  were 
applied,  and  the  result ingterms  were  combined,  the 
Patterson  flinction  was  written  as 

P-4/V  +  F^kil  cos(2'rrhu)cos(2Trkv)cos(2'n’lw) 

-  t'Skl-^’fci]  2^1  hu )  cos{2T['kv)sin(2Triw)  ^ . 
With  the  Patterson  function  in  this  form,  it  was  possible 
to  obtain  a  three-dimensional  Patterson  synthesis  using 
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the  BXFS  program.  Layers  of  the  Patterson  synthesis  upon 
which  important  peaks  were  located  were  mapped  for  closer 
study . 

5.  Interpreting  the  Patterson  Synthesis 
a.  The  Patterson  Projection.  P(u.w) 

The  map  of  the  zero-level  Patterson  projection 
showed  only  a  few  large,  well-defined  peaks,  each  of  these 
peaks  being  of  one  of  two  different  heights.  It  was  as¬ 
sumed  that  the  larger  of  the  two  kinds  of  peaks  represented 
a  bromine-bromine  interatomic  vector,  and  that  the  smaller 
peak  represented  a  sodium-bromine  interatomic  vector. 

Based  upon  these  assignments,  two  dimensional  coordinates 
for  the  bromine  and  sodium  atoms  were  found.  A  structure 
factor  calculation  showed  the  discrepancy  factor,  R,  to 
be  0.7;87,  (See  section  III  C3  for  a  discussion  of  this 
type  of  calculation  end  the  resulting  R-value.)  The  con¬ 
struction  of  the  zero-level  Patterson  projection  map  did 
serve  as  a  rough  guide  in  locating  bromine-bromine  peaks 
on  the  three-dimensional  Patterson  maps,  but  had  there 
been  computer  facilities  initially  available  for  the 
reedy  evaluation  of  the  three-dimensional  synthesis,  the 
projection  would  have  been  completely  unnecessary. 

b.  The  Three-Dimensional  Patterson  Synthesis 
The  advantage  of  employing  a  three-dimensional 
Patterson  synthesis  lies  in  that  the  resulting  peaks  (for 
a  simple  crystal)  usually  do  not  seriously  overlap. 

The  three-dimensional  Patterson  showed  three 


very  large  peaks,  in  addition  to  the  usual  origin  peak, 
which  were  located  at  u=25/60,  v»14/60,  w=35/60;  u-O/60, 
v»16/60,  w=30/60;  and  u»25/60,  v=30/60,  w=5/60, ,  It  was 
assumed  that,  since  bromine  has  by  far  the  largest  atomic 
number  of  the  elements  present  in  sodium  bromide  dihydrate, 
these  peaks  were  most  probably  representative  of  bromine- 
bromine  interatomic  vectors.  Tentative  coordinates  for 
the  bromine  atom  were  obtained  by  the  following  reasoning. 

The  general  equivalent  positions  of  the  space 
group  P2i/c  are  x,y,z;  x,4-y,^+z;  x,4*y,i-z;  and  x,y,2. 
Position  x,y,z  could  be  chosen  arbitrarily  as  a  starting 
point  for  formulating  predicted  interatomic  vectors.  The 
interatomic  vectors  l-2x,l-2y,l-2z ;  l-2x,^,i-22;  and 
0,i-2y,^  would  then  be  expected  to  appear  on  the  Patterson 
sections.  Now  if  the  peak  located  at  u=0/60,  v-16/60, 
w=30/60  corresponds  to  the  predicted  peak  at  0,i-2y,i,  then 
30/60-2y=16/60 ,  and  y-7/60.  In  like  manner,  if  the  peak 
located  at  u=25/60,  v=30/60,  w=5/60  corresponds  to  the  pre¬ 
dicted  peak  at  l-2x,|r,i-2z,  then  x  and  z  should  be  17.5/60 
and  12.5/60,  respectively.  Thus,  with  the  general  co¬ 
ordinates  of  the  bromine  atom  being  tentatively  established 
as  x-17.5/60,  y=7/60,  z=l2.5/60,  it  was  possible  to  predict 
that  the  third  bromine-bromine  vector  should  appear  at 
u-25/60,  v-14/60,  w*=35/60  (that  is,  at  l-2x,l-2y,l-2z) . 

The  third  of  the  peaks  did,  indeed,  appear  at  that  position. 
This  postulated  structure  was  tested  by  means  of  a  struc¬ 
ture  factor  calculation  as  explained  in  the  next  section. 
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C.  Locating  Atomic  Positions  In  the  Trial  Structure 

1,  The  Structure  Factor 

The  structure  factor  Fhklf  a  function  of  indices 
hkl,  may  be  written  in  general  as 

Fhkl  * fj  (hxj+kyj+lzj  ) .  The  scattering  factor,  f  J  , 

j  »( 

for  the  jth  atom  of  the  N  atoms  present  in  the  unit  cell, 
is  a  function  of  the  number  and  distribution  of  the  elec¬ 
trons  in  the  atom;  it  also  depends  upon  sin  ©/;(  ,  where  & 
is  the  Bragg  angle,  and  ^  is  the  wavelength  of  the  X-radi¬ 
ation  employed  in  the  diffraction  experiment.  The 
modulus  of  Fhkl -is  ”the  ratio  of  the  amplitude  of  the 
radiation  scattered  in  the  reflection  hkl  by  the  contents  of 
one  unit  cell  to  that  scattered  by  a  single  electron  under 
the  same  conditions”  (Lipson  and  Cochran,  195K,  p.  11). 

2.  The  Temperature  Factor 

When  discussing  a  crystal  structure,  it  must  be 
realized  that  the  atoms  involved,  far  from  remaining 
stationary,  undergo  thermally- induced  motion.  This  means 
simply  that  the  electrons  of  an  atom  are  ’’spread  out” 
over  a  larger  volume  in  space  than  they  would  occupy  if 
the  atom  were  at  rest.  The  scattering  factor  curve  for 
the  atom-in-thermal-motion  thus  decreases  more  rapidly, 
as  sin  ©/x  increases,  than  it  would  if  the  atoms  were  at 
rest.  An  approximate  correction  to  be  applied  to  the 
scattering  factor  is  the  Debye-V/aller  correction,  so  that 

(Buerger,  1960,  p.  231-2i52) ,  where  f 
is  the  modified  scattering  factor  for  a  given  atom,  is 


the  scattering  factor  for  that  stationary  atom,  and  B  is 
denoted  as  the  temperature  coefficient  or  as  the  isotropic 
temperature  factor  for  the  atom.  (B  is  generally  different 
for  each  of  the  different  kinds  of  atoms  in  a  structure.) 

Thermal  motion  may  be  considered  to  be  signifi¬ 
cantly  anisotropic  as  well  as  different  for  each  kind  of 
atom.  It  is  possible  to  convert  an  isotropic  temperature 
factor  into  its  corresponding  anisotropic  temperature 
factors  by  the  relation 

B(sin2e  ^12hk+  (?23kl+  (?3llh 

-  (Ba*V4)h2+(Bb*V4)k2+(Bc'^/4)l  +  ( Ba%*cos  r Vz ) 

+  iBb*c*coscX*/2)  +  (Bc*a*cos(?V2). 

By  equating  coefficients  of  h2,  k2,  12,  hk,  kl,  and  Ih,  one 
may  evaluate  the  ^’s  in  terms  of  B  and  the  usual  recip- 
rocal  cell  constants.  For  example,  (Ba*  /4) , 

(?22  -  (Bb*2/4)  , ...  These  relations  are  useful  in  that 
very  often,  for  purposes  of  refinement,  it  is  desired  to 
allow  the  anisotropic  factors  corresponding  to  a  known 
isotropic  factor  to  vary  independently,  as  in  a  least 
squares  refinement,  or  in  difference  map  refinement.  The 
variation  of  these  factors  will  depend,  of  course,  upon  the 
extent  to  vhich  the  atoms  undergo  anisotropic  thermal 
motion;  it  will  also  depend  upon  whether  or  not  certain 
undetected  or  uncorrected  systematic  errors  remain  in  the 
structure  factor  data. 

g.  The  Structure  Factor  Calculation  and  the  R-Value 
If,  for  a  particular  crystal  structure,  the  cal¬ 
culated  structure  factor  magnitudes  approximately  match  the 
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observed  structure  factor  magnitudes,  it  usually  indicates 
that  the  correct  model  for  that  crystal  structure  has  been 
proposed.  Therefore,  the  structure  factor  calculation  is 
useful  as  a  means  of  testing  v/hether  or  not  a  reasonable  ap¬ 
proximation  to  the  crystal  vStruct'jre  has  been  made. 

It  has  been  customary  to  regard  the  crystallo¬ 
graphic  discrepancy  factor,  R,  defined  as  R  =  SljFoMfcll 

^l^ol 

as  an  indication  of  the  correctness  of  a  structure.  Sub¬ 
scripts  ”0”  and  "c”  refer  to  observed  and  calculated  values 
of  |F|,  while  the  summation  is  taken  over  all  available  re¬ 
flections.  It  is  not  unusual  to  begin  refinement  of  a 
tentative  structure  wi.th  an  R-value  of  0.40,  and  to  work 
down  to  a  hi^ly-refined  structure  with  an  R-value  in  the 
neighborhood  of  0.05-0.10.  But  it  must  be  emphasized 
that  the  R-value  need  not  necessarily  be  accepted  at  its 
face  value.  Other  residuals  have  been  suggested  as 
measures  of  discrepancy,  among  them  being  Booth’s 
residual,  R  »  T  ,  shown  here  as  an  example. 

2,1^0! 

For  a  fuller  discussion  of  the  above,  see  Buerger’s  book 
(1960,  pp.  585-589).  In  this  study,  a  weighted  R-value 
R„  »  '^ilFol  -iFcll  is  calculated  in  addition  to  the 

!2.ui-'[7pV  2 

conventional  R-value.  The  weight,  w,  is  obtained  from 
the  equation  w  =  l/c^,  where  (T  is  the  estimated  standard 
deviation  in  1f|  .  Usually  in  judging  the  accuracy  of  a 
crystal  structure  model,  it  is  very  helpful  to  compare 
the  individual  calculated  structure  factors  with  their 
experimentally-obtained  values. 
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For  convenience  in  calculation,  the  form  of 
Fhki  (given  above) is  usually  modified  by  taking  advan¬ 
tage  of  the  space  group  symmetry  of  the  crystal,  as  is 
illustrated  below  for  P2i/c.  The  four  general  equivalent 
positions  of  this  space  group  are  x,y,2;  x,i-y,^'-*-z;  x,i+y,|-z 
5c, y, 2.  The  structure  factor,  may  then  be  written  as 

Fhkl  ■  X  1. 

J  *  \  ^*1 

and  oC  refers  to  one  cd'  the  number  of  atoms  of  kind  j.  The 
summation  over  may  be  evaluated  in  the  following  manner: 

*^Aj,oc,hkl  =  e2‘"i(hxj.kyj+lzj)  ^  ^-STT  i(hXj+kyj-vlZj ) 

+  gf  Tri(hxj+k(i-yj)  +  l(i+zj) ) 

^  -2Tri(hxj+k(i-yj )+l{i+zj) . 

Since  the  general  relation,  ©six  *cos(ax)+isin(ax)  holds, 

the  preceding  expression  may  be  simplified  to  yield 

tAj.or.hkl  =  2cos(2TT^hxj+kyj+lzjl) 

1 

+2cos(2TtChXj-kyj+lZj+(k-vl)/2 ) ) , 

Appl? ''■=tion  of  the  identity  cos®C +008^  =2cos(^(o(+^) ) 
cos  (i(Qc' -^  ) )  produces  Aj  ,c5f  ,hkl=  4|cosi(4'Tt  (hx j+lzj ) 

+2Tt  (k+l)/2)cos^(4Tf  kyj-2Tl’ (k+l/2)|  ,  which  may  be  simplified 
to  4cos2'n' (hxj+lzj+(k+l)/4)cos2Tr(kyj-(k+l)/4) .  Thus  the 
structure  factor  is  Fhkl  “  ^  ^  f  jC0s2Tt' {hxj-wlzj+(k+l)/4) 
cos2Tt  (kyj-(k'^l)/4) .  A  further  simplification  may  be  ob¬ 
tained  by  splitting  the  summation  into  two  portions,  de¬ 
pending  upon  whether  (k+l)  is  even  or  odd,  then  applying 
simple  trigonometric  identities. 
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The  results  are 

^hkl  ®  4  X  ,  f  i  |lCos{2‘TT  hxj)cos(2TT  kyj  )cos(2'TTl2j  ) 
-sin(2^hxj)cos(2TTkyj)sin{2Ti  Izj)  j  ,  and 

FHkl  *  -4  f  j  [sin{2Tt  hxj)sin(2Tt' kyj  )cos  (2TT  Izj ) 

J - 

+  COS  (2TT  hxj  )sin{2Tr  kyj  )sin(  211  Izj  )  j  . 

In  these  mathematical  forms,  the  values  of  the  structure  factors 
may  be  efficiently  calculated  using  standard  tables  of 
sin(2TThx)  and  cos{2Trhx)  (International  Tables,  1959,  Sec. 8). 

In  practice,  however,  all  structure  factor  calcu¬ 
lations  were  performed  by  computer  means.  computer  prog¬ 

rams  were  utilized — the  Fittsburgh  Structure  Factor  Program 
(Shiono,  1962),  a  means  to  calculate  rapidly  structure  factors 
with  isotropic  temperature  factors  only,  and  the  BXLS  Pro¬ 
gram  (Carpenter,  1963),  which  facilitates  the  calculation  of 
structure  factors  with  anisotropic  temperature  parameters.  The 
BXLS  program  output  also  provides  values  of  (Fq-Fq)  for  all 
hkl  reflections,  which  may  thai  be  used  as  Fourier  coeffi¬ 
cients  for  a  difference  Fourier  synthesis  if  so  desired. 

Tables  of  atomic  scattering  factors,  obtained  by  quantum- 
mechanical  calculations  from  self-consistent  or  variational 
wave  functions  are  readily  available  for  the  sodium  cation,, 
bromine  anion  and  oxygen  atom  (International  Tables,  1962, 
pp.  202-207). 

A  structure  factor  calculation  was  made  using 
only  the  tentative  coordinates  of  the  bromine  atom;  an  iso¬ 
tropic  temperature  factor  parameter  for  bromine  was 
assumed  to  be  2.0A‘^.  The  R-val\ie  of  the  calculation  was 

0.258,  indicating  that  a  reasonable  model  had  indeed 
been  proposed,  also  that  a  good  deal  of  the  X-ray 
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scattering  in  the  unit  cell  of  sodium  bromide  dihydrate  is 
due  to  the  bromine  ion  only.  From  the  calculation,  the 
tentative  assignment  of  phases  (positive  or  negative  signs 
in  this  case)  to  the  observed  structure  factors  could  be 
made.  Knowledge  of  both  the  magnitudes  and  the  phases  of 
the  observed  structure  factors  made  possible  the  construc¬ 
tion  of  a  three-dimensional  electron  density  map. 

4.  The  Electron  Density  Calculation 
As  mentioned  earlier,  the  general  expression  for  the 
electron  density  evaluated  at  point  x,y,z  in  the  unit  cell 
is  p(x,y.i).  (Lipson  and 

Cochran,  1953,  p.l2),  the  summation  being  performed  over 
indices  hkl.  Since  the  crystal  under  investigation  is 
centosymraetric ,  the  form  P(x,y,z)  =  1/V  /  /  /  Fhkl 

—  M 

cos2Tf{hx+ky4-lz )  is  easily  obtained,  as  in  the  treatment  of 
the  three-dimensional  Patterson  function  (see  III  Bl).  The 
sum  was  split  into  two  portions,  one  for  the  case  of  in¬ 
dices  k+1  =  2n,  and  another  for  the  case  of  indices  k+1  = 
2n+l.  Use  of  the  identities  F^kl*  ^hTcli  FRkl  “  ^‘‘hkl  fpr 
k+1  =  2n,  and  F^kl  =  -Fhkl,  Fj^ki  =  -Fhkl  for  k+1  -  2n+l, 
coupled  with  the  application  of  standard  trigonometric 
identities,  led  to  the  following  form  of  the  electron 
density  function  suitable  for  computer  calculation. 
f(x,y,z)-4/v21  2.'2^kl+"FRkl3  c  0  s  ( 2TThx )  s  i  n  ( BTT’ky )  s  i  n  ( 2Ttl  z ) 

^  +  lj^ikl-Fhk3sin(2TThx)  sin(2rrky)cos(2Triz)J' 

-4/V  /  2l»\L^hkl+Fhkl3cos(2rrhx)sin(2rrky)sin(2n'lz) 

+  [Fhkl-FHkl3  5in(2nhx)sin(2n’ky)cos(2rrlz)^  , 
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5 .  Obtaining  Atomic  Coordinates . 

Four  separate  peaks  appeared  on  the  first  three- 
dimensional  electron  density  map.  By  examination  of  the 
relative  peak  heights,  it  was  possible  to  label  the  peaks 
unambiguously  as  bromine,  sodium  and  the  two  oxygens.  It 
was  found  convenient  to  locate  peak  canters  by  an  ap¬ 
proximation  method  (Carpenter  and  Donahue,  1950).  The 
coordinates  of  these  peak  eerters  were  taken  as  the  ten¬ 
tative  coordinates  of  the  four  atoms.  After  a  second 
structure  factor  calculation  was  made,  using  these  co¬ 
ordinates  together  with  an  assumed  isotropic  temperature 
°2 

factor  of  2.0A  for  all  the  atoms,  the  R-value  was  a  much- 
improved  0.127.  This  second  calculation  provided  an 
improved  set  of  structure  factor  signs  for  a  second 
electron  density  calculation  which,  in  turn,  provided  im¬ 
proved  atomic  coordinates  for  a  third  stinicture  factor 
calculation.  After  the  results  of  the  third  calculation 
were  examined,  it  was  apparent  that  there  were  no  further 
changes  to  he  made  in  the  assignments  of  phases  to  the 
structure  factors.  (The  above  procedure  is  often  referred 
to  as  the  method  of  refinement  by  Fourier  cycles 
(Buerger,  1960,  pp.  590-593),)  A  few  more  factor  cal¬ 
culations  led  to  the  adjustment  of  the  overall  isotropic 
temperature  factor,  from  its  initial  value  of  2.0A^  to 

o2 

the  value  2.5A  .  As  indicated  by  the  R-value  of  0.103, 
there  was,  at  this  point,  fair  agreement  between  calcu¬ 
lated  and  observed  structure  factors.  In  order  to 
refine  the  structure  further,  it  was  decided  to  employ  the 
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three-dimensional  difference  synthesis. 

D.  Refining  the  Trial  Structure 

1.  The  Difference  Electron  Density  Calculation 
The  general  expression  for  the  electron  density 
using  the  experimentally-derived  structure  factor  values 
as  Fourier  coefficients  may  be  written  as 
fo  =  l/V^, i(hx+ky->-l^z)  ^  y,here  subscript  ”o” 
indicated  "observed”.  If  structure  factors  based  upon 
a  proposed  model  oi'  the  crystal  structure  are  used  as 
Fourier  coefficients,  the  electron  density  may  also  be 
expressed  as  fg  =  l/V  ^ Fge"^^  ,  where 

subscript  ”c"  indicates  "calculated".  If  Pc  is  sub¬ 
tracted  from  Pq,  the  expression 

(Po-p.)-iA2rz  (Fo-Fc  results; 

it  is  a  general  expression  of  the  difference  electron  den¬ 
sity  function,  utilizing  (Fq-Fq)’s  as  Fourier  coefficients. 
The  form  of  (  P o”  suitable  for  computer  calculation 

may  be  obtained  by  employing  exactly  the  same  procedure 
as  was  used  in  dealing  with  the  electron  density  function. 

If  the  trial  structure  exactly  matches  the 
real  structure,  the  difference  map  should  be  absolutely 
featureless,  but,  in  reality,  random  fluctuations  will 
appear  owing  to  observational  errors.  If,  however,  the 
trial  structure  differs  in  some  systematic  way  from  the 
real  structure,  the  difference  map  will  often  exhibit 
that  difference,  and,  in  so  doing,  will  indicate  the 
direction  in  which  improvements  of  atomic  parameters 
should  be  carried  out  in  order  to  produce  a  more  nearly 
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featureless  map. 

2.  Cerrectien  of  Caordlnate  Errors  by 
Difference  Maps 

If  a  postulated  atomic  location  is  only 

moderately  in  error,  the  difference  map  will  reveal  it  to 

be  situated  upon  a  steep  gradient.  There  is  a  convenient 

semi-empirical  formula  used  to  evaluate  the  magnitude  of 

the  correction.  It  is  the  relation 

€  =  d(Po  -  Pc)  /  2  (Buerger,  1960,  pp.  604-606), 

where  g  is  the  coordinate  correction,  r  is  the  distance 

from  the  atomic  center,  d(  Pp-  Pg)  is  the  slope  of  (Pq-  P#) 

dr 

in  the  immediate  neighborhood  of  the  center,  and  Pg(0) 
is  the  maximum  electron  density  given  approximately  by 
Po(0)  =  Z(p/'n'  Z  denoting  the  atomic  number  and  p 

being  about  5.0,  The  fractional  coordinate  corrections 
carried  out  for  the  bromine,  sodium  and  two  oxygen  atoms 
were  all  quite  small,  none  exceeding  0,002, 

3,  Ad.iustment  of  Anisotropic  Temperature  Factors  by 

Difference  Maps 

On  the  difference  map,  anisotropy  in  the  thermal 
motion  of  an  atom  will  generally  be  evidenced  by  the  ap¬ 
pearance  of  a  saddle-shaped  region  about  the  atomic 
center.  (As  previously  mentioned,  this  appearance  might 
also  be  traced  to  the  presence  of  systematic  errors  in 
the  structure  factor  data.)  Some  relatively  large  adjust¬ 
ments  had  to  be  made  in  the  temperature  parameters  of 
bromine,  sodium  and  the  two  oxygens.  After  adjustments 
based  upon  four  difference  maps  had  been  carried  out,  it 


was  jvidged  that  the  expense  of  computer  operation  and 
the  time  involved  in  the  determination  of  smaller  and 
smaller  parameter  adjustments  seemed  not  to  warrant 
further  refinement  by  this  method.  A  good  general 
discussion  of  the  difference  synthesis  is  given  by 
Buerger  (1960,  pp.  603-609)  and  Lipson  and  Cochran  (1953, 
pp.  208-306). 


For  the  final  refinement  of  the  crystal 
structure,  a  least  squares  method  seemed  most  suitable. 

4.  Least  Squares  Refinement 

The  BXLS  computer  program,  which  carries  out 
refinement  by  least  squares,  minimizes  essentially  the 

vr  1  KlFdi  -  lFc\l  ^  , 

function  wIKIFqII^  ~  t )•  The 

minimization  of  R^i  led  to  changes  in  all  of  the  atomic 
parameters.  After  two  cycles  of  least  squares  were  com¬ 
pleted,  a  difference  map  was  constructed.  The  map  showed 
marked  improvement  in  that  it  was  much  more  nearly 
featureless  than  previous  ones.  After  the  third  cycle  of 
least  squares,  however,  it  was  evident  that  the  parameter 
values  had  converged,  that  is,  the  changes  in  the  para¬ 
meters  were  insignificantly  small.  At  that  point  it  was 
judged  that  the  refinement  was  completed. 

A  total  of  506  reflections  were  employed 
throughout  the  course  of  refinements;  no  "unobserved"  re¬ 
flections  were  utilized.  (an  unobserved  reflection  is 
one  which  has  a  corresponding  intensity  value  less  than  a 
certain  observable  minimum.)  The  comparisons  between 
the  observed  and  ealculeted  structure  factors,  for  all 


the  available  hkl  reflections,  are  listed  in  the  Appendix 
of  this  thesis  . 

The  effects  of  the  difference  map  and  least- 
squares  refinements  upon  the  atomic  parameters  and  other 
structural  quantities  which  resulted  from  Fourier  refine¬ 
ment,  are  diovrn  in  Table  II.  It  should  be  remembered 
that  the  difference  map  refinement  was  not  carried  as  far 
as  it  should  have  been. 
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TABLE  II 

Effects  of  Difference  Map  and  Least-Squares  Refinement  upon 
Atomic  Parameters  and  Other  Structural  Quantities 

Bromine 


Parameter 

Fourier 

Refinement 

Diff.  Map 
Refinement 

Least  Sq. 
Refinement 

X 

0.2962 

0.2969 

0.2964 

y 

0.1194 

0.1195 

0.1191 

z 

0.2117 

0.2123 

0.2119 

^  11 

0.01718 

0.02406 

0.02404 

^22 

0.00572 

0.00800 

0.00841 

^33 

0.01615 

0.01939 

0.01861 

^12 

0.0 

0.0 

0.0 

^23 

0.0 

0.0 

0.0 

^31 

0.01322 

0,01851 

0.01875 

X 

0.0150 

Sodium 

0.0146 

0.0146 

y 

0.1680 

0,1684 

0.1693 

z 

0.4567 

0,4568 

0.-4567 

(?11 

0.01718 

0.02406 

0.02525 

^22 

0.00572 

0.00800 

0.00772 

^33 

0.01615 

0,01939 

0.01798 

^12 

0.0 

0.0 

0.0 

^23 

0.0 

0.0 

0.0 

^31 

0,01322 

0.01851 

0.01263 

TABLE  II  (continuad) 
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Oxygen 

(I) 

Parameter 

Fourier 

Refinement 

Diff.  Map 
Refinement 

Least  Sq. 
Refinement 

X 

0.7940 

0.7940 

0.7882 

y 

0.3110 

0.3118 

0.3132 

7, 

0.1780 

0.1774 

0.1762 

(?11 

0.01718 

0.02268 

0.02052 

^22 

0.00572 

0,00755 

0.00831 

^33 

0.01615 

0.02132 

0.02213 

^12 

0.0 

0.0 

0.0 

(^23 

0.0 

0.0 

0.0 

^31 

0.01322 

0.01746 

0.00911 

Oxygen 

(ID 

X 

0.2129 

0.2139 

0.2130 

y 

0.4929 

0,4914 

0.4900 

z 

0.2218 

0,2203 

0^2199 

0.01718 

0.02268 

0.02063 

(^22 

0.00572 

0.00755 

0.00833 

^33 

0.01615 

0.02132 

0.02292 

(^12 

0.0 

0.0 

0.0 

^23 

0.0 

0.0 

0..0 

P3I 

0.01322 

0.01746 

0.01608 

TABLE  II  (continued) 


Quantity 

Other 

Fourier 

Refinement 

Quantities 

Diff.  Map 
Refinement 

Least  Sq, 
Refinement 

K( scale 

4.40469 

3.91000 

3.89510 

factor ) 

R- value 

0.103 

0.065 

0.0615 

Rvf- value 

0.133 

0.077 

0.0707 

Note:  The  anisotropic  parameters  ^23 

required  to  be  equal  to  zero  since,  for  the  mono 
clinic  crystal,  cos  of*  =  cost’*'  =0  (see  III  C  2) 
Therefore,  these  parameters  were  held  invariant 
throughout  the  structure  refinement. 
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IV.  Description  and  Discussion  of  the  Structure 

1 

A.  General  Structural  Features 
« 

Figure  1  is  a  projection,  on  the  a-£  plane,  of 
the  crystal  structure  of  sodium  bromide  dihydrate.  The 
numbers  in  parentheses  are  fractional  coordinates  which 
indicate  distances  along  the  b  axis.  The  structure  may 
be  described  as  composed  of  separate  layers  parallel  to 
the  h-c  plane.  Each  of  these  layers  may  be  subdivided 
into  three  sheets  of  atoms  containing,  respectively, 
bromine  end  oxygen,  sodium,  and  bromine  and  oxygen.  As 
will  be  more  fully  discussed  later,  the  separate  layers 
appear  to  be  held  together  largely  by  hydrogen  bonds, 
thus  accounting  for  the  ready  cleavage  of  this  hydrate  in 
the  b-c  plane  (Van  Meersche  et  al,  1962). 

Alternately,  one  may  regard  the  layers  as  being 
constructed  from  a  series  of  octshedra,  with  the  sodium 
cation  at  the  center  of  each  octahedron.  The  sodium 
cation  is  coordinated  to  four  oxygen  atoms  and  to  two 
bromine  anions,  but  each  bromine  anion  and  each  oxygen 
atom  is  also  coordinated  to  an  adjoining  sodium  cation. 
Thus,  since  each  electxx)negative  atom  is  shared  between 
two  adjoining  octahedra,  the  cohesion  within  the  layers 
parallel  to  the  b-c  plane  is  explained  (see  Figure  2A). 

*  The  figures  are  on  pages  54  56. 


44. 


B.  Discussion  of  Interatomic  Distances 

Based  upon  the  coordinates  of  the  sodium, 
bromine  and  tw  oxygen  atoms  \ihich  resulted  from  least- 
squares  refinement,  and  the  unit  cell  dimensions,  it  was 
possible,  with  the  aid  of  an  ’’Interatomic  Distances”  com¬ 
puter  program  (Carpenter  and  Hall,  1962)  to  obtain  all  the 

V  •  ■ 

relevant  interatomic  distances  in  the  sodium  bromide  dihydrate 
crystal.  These  distances  are  listed  below  with  their  stand¬ 


ard  deviations  in  parentheses.  The  standard  deviations  were 
calculated  from  the  equation  for  variance,  <r^ 
where  D12  refers  to  the  distance  between  atom  1  and  atom  2, 
is  the  partial  derivative  with  respect  to  the  ith 
coordinate.  (The  variance,  is  the  square  of  the 

standard  deviation  obtained  from  BXLS  output;  it  is  an  esti¬ 
mate  of  internal  data  consistency  only.) 

The  distances  of  2.965A  (0.005A)  and  2.98?A(0,006A) 


between  the  sodium  and  bromine  ions  are  consistent  with  the 
sum  of  their  crystal  radii,  2.90X  (Pauling,  1960,  p.  514). 


The  two  distances  are  almost  identical  with  that  observed 
in  the  sodium  bromide  crystal,  2.98^.  The  sodium  ion- 
oxygen  atom  distances  of  2.40lX(0.009X) ,  2.410A(0.009X) , 
2.418X(0.009X) ,  and  2.495X(0.1lX)  are  also  in  good  accord 
with  the  predicted  value  of  2.35X  based  upon  Pauling's 
crystal  radii  (see  Figure  2A). 


The  shortest  oxygen-bromine  ion  distances  are 
3.355A(0.008A) ,  3.356X(0.009A),  and  3.383X(0.010A) .  The 

O 

sum  of  the  van  der  Waals  radii  of  these  atoms  is  3,35A 


(Pauling,  1960,  p.  260).  The  next  shortest  bromine-oxygen 
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interatomic  distance  is  3.575X{0.007X) .  An  important 
oxygen-bromine-oxygen  angle  is  103.35“{0.27®)  (see  Figure  1), 
The  following  interatomic  distances  and  angles 
refer  only  to  atoms  which  belong  to  one  octahedron. 

O  O 

Bromine-oxygen  distances  vary  from  3.806A(0.009A )  to 

0  0  e 

3.920A{0.009A) ,  oxygen-oxygen  lengths  range  from  3.185A 
(O.OOqA)  to  3.637X(0.013X) ,  and  the  two  bromine  ions  are 
4.356^(0.002^)  apart.  The  angle  formed  by  the  sodium 
cation  and  the  two  bromine  anions  is  94,22f(0.i7" ) ,  while 
the  oxygen- sodium- oxygen  and  oxygen-sodiura-bromine  angles 
vary  from  82.92* (0.42® )  to  95.50"(0.57®) (see  Figure  2B) . 

C,  Comparison  Between  This  Structure  Determination 
and  Van  Meerssche^s 

The  structure  of  sod>;ium  bromide  dihydrate  just 
described  is  almost  exactly  that  worked  out  by  Van  Meerssche 
and  co-workers  (1962).  This  development  presented  a  rare 
opportunity  for  a  comparison  of  two  completely  independent, 
yet  contemporary,  structure  determinations  of  the  same  sub¬ 
stance.  There  is,  in  general,  good  agreement  between  the 
atcsnic  coordinates  determined  here  and  those  reported  by 
Van  Meerssche  and  co-workers  at  the  University  de  Louvain 
in  Belgium.  But  in  order  to  arrive  at  a  more  detailed 
comparison  of  the  two  structure  determinations,  it  was 
decided  to  subject  Van  Meerssche’ s  data  to  the  BXLS 
least-squares  program  previously  mentioned.  In  this 
treatment,  the  weights  assigned  to  Van  Meerssche ’s  struc¬ 
ture  factor  data  were  bared  upon  those  employed  in  this 
study.  In  addition,  his  isotropic  temperature  factors 
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were  converted  to  anisotropic  temperature  factors,  which 
were  then  allowed  to  vary  independently  along  with  the 
coordinates  and  scale  factor.  In  essence,  the  final  re¬ 
finements  of  the  two  sets  of  data  resulting  from  the  two 
separate  studies  were  identical.  Table  III  is  presented 
fcr  purposes  of  comparison. 

It  may  be  readily  seen  that  by  allowing  each  of 
the  anisotropic  temperature  factors  to  vary,  slight 
changes  appear  in  the  coordinates  originally  reported  by 
Van  Meerssche.  Indeed,  these  changes,  on  the  whole, 
lead  to  less  agreement  between  the  two  sets  of  coordinates. 

In  most  cases,  however,  the  altered  coordinates  do  not 
differ  from  those  obtained  in  this  study  by  more  than  2(r  , 
where  0*  is  an  estimated  standard  deviation  obtained  from 
the  BXLS  program  output.  To  be  more  specific:  Of  the 
twelve  coordinates,  six  agree  within  the  range  of  0  -T , 
four  within  d"  -20"  ,  and  two  within  2  (T  -  36“  , 

No  useful  comparison,  however,  can  be  made 
between  the  two  sets  of  anisotropic  temperature  factors, 
owing  to  the  different  absorption  errors  in  both  sets  of 
structure  factor  data.  This  point  is  more  fully  discussed 
in  the  following  paragraphs. 


TABLE  III 
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Structural  Parameters  from  Brown  University  and  from  the 
Universite  de  Louvain,  after  BXLS  Refinement 

Bromine 

Note:  The  figures  in  pai'entheses  are  standard  deviations 
obtained  from  BXLS  output. 

University  de  Louvain 


Parameter 

Brown  Universitv 

retorted 

after  BXLS 

X 

0.2964  (0.00017) 

0.2965 

0.2965  (0.00020) 

y 

0.1191  (0.00013) 

0.1191 

0.1192  (0.00014) 

2 

0.2119  (0.00018) 

0.2120 

0.2120  (0.00021) 

(?11 

0.02404(0.00024) 

0.02197 

(B-3.19) 

0.01995(0.00032) 

^22 

0.00841(0.00012) 

0.00740 

0.00704(0.00015) 

^33 

0.01861(0.00022) 

0.02063 

0.02020(0.00033) 

^12 

0.00000(0.00000) 

0.00000 

0.00000(0.00000) 

0.00000(0.00000) 

0.00000 

0.00000(0.00000) 

^31 

0.01875(0.00034)  0.01698 

Sodium 

0.01539(0.00049) 

X 

0.0146  (0.00067) 

0.0155 

0.0159  (0.00075) 

y 

0.1693  (0.00049) 

0.1679 

0.1683  (0.00050) 

z 

0.4566(0.00069) 

0.4576 

0.4579  (0.00075) 

<?ll 

0.02525(0.00117) 

0.02169 

(B-3.15) 

0.02198(0.00116) 

^22 

0.00772(0.00050) 

0.00731 

0.00465(0.00045) 

^33 

0.01798(0.00108) 

0.02037 

0.01748(0.00113) 

(?12 

0.00000(0.00000) 

0.00000 

0.00000(0.00000) 

^23 

0.00000(0.00000) 

0.00000 

0.00000(0.00000) 

^31 

0.01263(0.00180) 

0.01676 

0.01610(0.00173) 

TABLE  lil  (continued) 


Oxygen  (I) 


Parameter 

Brown  University 

University  de  Louvain 
reported  after  BXLS 

X 

0.7882  (0.00122) 

0.7907 

0.7915  (0.00131) 

y 

0.3132  (0.00091) 

0.3132 

0.3138  (0.00102) 

z 

0.1762  (0.00129) 

0.1761 

0.1770  (0.00149) 

^11 

0.02052(0.00184) 

0.02472 

0.01924(0.00209) 

^22 

0.00831(0.00095) 

(B-3.59) 

0.00833 

0.00580(0.00091) 

^33 

0.02213(0.00211) 

0.02321 

0.02367(0.00233) 

^12 

0.00000(0.00000) 

0.00000 

0.00000(0.00000) 

^23 

0.00000(0.00000) 

0.00000 

0.00000(0.000  00) 

^31 

0.00911(0.00302) 

0.01911 

0.01497(0.00328) 

X 

Oxygen 

0.2130  (0.00112) 

(II) 

0.2162 

0.2154  (0.00132) 

y 

0.4900  (0.00084) 

0.4922 

0.4919  (0.00096) 

z: 

0.2199  (0.00130) 

0.2195 

0.2189  (0.00143) 

^11 

0.02063(0.00175) 

0.02370 

0.02100(0.00210) 

^22 

0.00833(0.00099) 

(B=3.47) 

0.00805 

0.00651(0.00094) 

^33 

0.02292(0.00198) 

0.02244 

0.01974(0.00208) 

^12 

0.00000(0.00000) 

0.00000 

0.00000(0.00000) 

(^23 

0.00000(0.00000) 

0.00000 

0.00000(0.00000) 

^31 

0.01608(0.00284) 

0.01847 

0.01662(0.00320) 

49 


TABLE  III  (continued) 

Quantity 

Brown  University 

Universitl  de  Louvain 
reported  after  BXLS 

R( scale  factor) 

3.89510 

0.80 

0.88290 

R(obs.only) 

0.0615 

0.0999 

Rjobs.only) 

0.0707 

0.1188 

R(all  refl.) 

0.1236 

0.1095 

R(all  refl. ) 

w 

(0.114) 

0.1486 

0.1239 

Shortest  0-Br 

3.355A 

3.32A 

3.325A 

Distances 

3.356X 

3.37A 

3.372A 

3.383X 

3.36A 

3.3621 

Notes:  ^12  and  (^25  ®nst  be  zero,  (see  Note  at  end  of 

Table  II. 

The  "reported"  R-yalues  of  the  UniversitI  de 
Louvain  are  jshose  which  resulted  from  a 
structure  factor  calculation  using  the 
"reported”  coordinates  and  B-values,  together 
with  the  BXLS  program.  The  R-value  in 
parentheses  is  taken  directly  from  Van 
Meerssche’s  paper  (1962). 


When  the  two  sets  of  coordinates  were  used  to 
calculate  structure  factors,  a  comparison  between  the  ob¬ 
served  and  calculated  structure  factors  brought  forth  two 
interesting  observations.  Firstly,  Van  Meerssche’s  group 
has  taken  pains  in  obtaining  intensity  data  for  very  weak 
reflections;  secondly,  in  spite  of  these  apparently  pains¬ 
taking  intensity  measurements,  with  an  integrating  Weis- 
senberg  goniometer,  the  agreement  between  observed  and  cal¬ 
culated  structure  factors  was  not  quite  as  good  (for  the 
same  reflections)  as  it  was  in  this  study  (Van  Meersche  et 
al,  1962;  also  see  Appendix). 

A  possible  clue  to  the  explanation  of  the  latter 
observation  may  be  found  in  the  fact  that  no  absorption  cor 
rections  were  applied  either  in  this  work  or  in  that  per¬ 
formed  by  Van  Meerssche’s  group.  The  absorption  correction 
would  be  more  important  for  data  obtained  from  the  rectangu 
lar  crystal  and  CuKo<.  radiation  employed  by  Van  Meerssche 
than  for  the  cylindrical  crystal  and  MoKo<  radiation  used 
in  this  thesis  work.  This  absorption  effect  is  notice¬ 
able  in  that  Van  Meersche ’s  temperature  factors  tend  to  be 
lower  than  those  obtained  here.  It  appears  that  the 
higher  precision  in  intensity  measurements  gained  by  Van 
Meerssche  and  co-workers  by  employing  an  integrating  Weis- 
senberg  camera  was  offset  by  their  omitting  the  absorption 
correction. 

The  structure  of  sodium  bromide  dihydrate 
proves  to  be  isomorphous  also  with  that  of  the  dihydrate 
of  sodium  cyanide  (Le  Blhan,  1958). 
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D.  Hydrogen  Bonding  in  Sodlmn  Bromide  Dlhvdrate 
Only  three  oxygen-bromine  interatomic  dis- 

o  o  o 

tances,  I5.355A,  b.SSbA  and  3,383A,  are  short  enough  to 
be  reasonably  indicative  cf  hydrogen  bonds.  Apparently, 
only  three  of  the  four  hydrogen  atoms  present  in  the  unit 
of  sodium  bromide  dihydrate  are  involved  in  hydrogen 
bonding.  This  finding  is  supported  by  Hornig's  infrared 
investigations  (see  Chap.  I)  and  by  Van  Meerssche's  nuclear 
magnetic  resonance  experiments.  The  latter  indicate  that 
nearly-rectilinear  hydrogen  bonds  exist  between  those  oxygen 

O  O 

and  bromine  atoms  separated  by  distances  of  3.32A,  3.37A 
and  3,36^  (Van  Meerssche ,  1962),  The  six  0~Br“  distances 
just  mentioned  are  very  similar  to  the  hydrogen  bonded 
0-Br"  distances  of  3,30X  and  3.35X  inD(-)  isoleucine 
HBr‘H20,  and  to  the  distance  of  3,38A  in  11-amino-undecanoic 
acid  HBr*i^H20  (quoted  by  Pimental  and  McClellan,  1960,  p. 
290)  . 

In  this  study,  it  was  determined  to  investigate 
the  possibility  of  locating  the  hydrogen  atoms  directly  by 
examination  of  difference  maps.  It  was  possible  to  detect 
three  peaks  on  the  difference  map  which  are  located  at  about 
the  positions  where  one  would  expect  to  find  the  hydrogen 
atoms  involved  in  hydrogen  bonding.  The  proposed  loca¬ 
tions  for  the  three  hydrogen  atoms  are  x  »  0.317,  y  =  0.567, 
z  =  0.250;  X  =  0.283,  y  »  0.422,  z  =  0.317;  and  x  =  0.733, 
y  -  0.395,  z  »  0,200,  The  oxygen-to-peak  distances  were 
calculated  to  be  0.96X,  1.02&,  and  0.9VA. .  The  H-O-H 
angle  was  then  calculated  to  be  110,3“  (see  Figure  1). 
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These  numbers  are  in  fair  agreement  with  the  known  para¬ 
meters  of  the  deuterium  oxide  molecule.  From  neutron 

0 

diffraction  work,  the  0-D  distance  is  listed  as  I.OIA,  and 
the  D-O-D  angle  is  givai  as  109.5“  (Pauling,  1960,  p.468). 

Unfortunately,  there  are  many  peaks  appearing  on 
the  difference  map  which  are  as  large  as^  or  slightly 
larger  than, the  three  peaks  just  mentioned.  Therefore,  al¬ 
though  the  three  peaks  in  all  probability  indicate  hydrogen 
atom  positions,  the  peaks  cannot  be  unambiguously  idenfied 
as  such. 

There  was  no  peak  discernible  on  the  difference 
map  which  could  be  designated  as  the  position  of  the  fourth 
hydrogen  atom.  Its  location  may  even  be  a  disordered  one: 
even  though  the  H-O-H  angle  must  be  maintained,  there  are 
many  possible  locations  for  the  fourth,  non-hydrogen-bonded 
hydrogen  atom. 

The  extent  of  hydrogen  bonding  in  sodium  bromide 
dihydrate  is,  oddly  enough,  much  greater  than  that  reported 
to  exist  in  the  dihydrate  of  sodium  cyanide  (,Le  Bihan, 

1958).  The  nitrogen  atom  typically  shows  a  markedly  greater 
tendency  towards  hydrogen  bond  formation  than  does  the 
bromine  atom.  Hence,  one  would  expect  that  at  least  three 
of  the  four  hydrogen  atoms  in  the  cyanide  (structurally 
isomorphous  with  the  bromide)  would  be  involved  in  hydrogen 
bonds.  But  Le  Bihan  found  only  one  nitrogen-oxygen  inter¬ 
atomic  distance  (2.8ll)  short  enough  to  be  considered  for 
hydrogen  bonding.  It  should  be  noted, how ever,  that  the 
R-value  for  the  cyanide  structure  is  0.19.  It  is 
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therefore  difficult  to  assess  the  reliability  of  Le  Bihan’s 
reported  interatomic  distances  for  the  only  partially- 
refined  structure.  For  example,  there  is  a  nitrogen-oxygen 

O 

distance  of  3.18A  in  the  structure  which,  after  further 
structure  refinement,  mi^t  possibly  be  modified  to  a 
value  indicative  of  at  least  a  weak  hydrogen  bond. 

In  summary,  this  thesis  work  has  shown  that 
a  well-refined  crystal  structure,  even  one  composed  of  both 
light  and  heavy  atoms,  can  provide  interesting  information 
regarding  the  distribution  of  hydrogen  bonds  in  the  struc- 
tiire,  and  the  probable  locations  of  the  hydrogen  atoms.  It 
was  gratifying  to  note  the  agreement  between  this  crystal 
structure  analysis,  the  infrared  work  done  by  Hornig  (1961), 
and  the  structure  determination  carried  out  by  Van 
Meerssche  and  co-workers  (1962). 


Figure  2 A 

Projection  of  the  NaBr'2H20 
structure  along  the  a  axis 
(only  2  octahedra  shovm) 
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Appendix 

Structure  Factors 
Obtained  In  This 
Thesis  Work 

(after  BXLS  refinement) 


Structure  Factors 
Qd.Ten  by  Van  Mserssche 
and  Co-workers  (1962) 
(after  BXLS  reflnanent) 


Notes  t  All  F  's  prefixed  by  a  minus  sign  are  "unobserred*. 
All  suructare  factors  are  anltlpUed  Iqr  10  for 
convenience. 


hkl 

^0 

'•c 

^c 

100 

184 

135 

159 

135 

200 

1057 

-  982 

865 

-  991 

300 

666 

659 

653 

687 

400 

528 

529 

600 

548 

500 

213 

-  216 

274 

-  256 

600 

104 

101 

124 

118 

700 

150 

154 

194 

197 

002 

1048 

-1043 

945 

-1036 

102 

104 

99 

115 

96 

T02 

767 

879 

962 

881 

202 

1202 

1212 

1210 

1234 

202 

958 

981 

848 

985 

2P2 

-  89 

0 

-  18 

.  4 

302 

434 

-  468 

397 

-  472 

402 

303 

-  311 

362 

-  336 

1102 

104 

100 

88 

115 

i02 

^^rZ 

275 

344 

311 

502 

ikZ5 

436 

477 

477 

602 

92 

128 

141 

141 

Zoz 

206 

-  200 

247 

-  225 

Zoz 

702 

.  142 

100 

-  69 

-  80 

106 

-  104 

802 

185 

185 

004 

550 

532 

486 

536 

104 

432 

437 

406 

445 

To4 

604 

-  584 

433 

-  559 

204 

394 

-  396 

406 

-  404 

.204 

-  32 

4 

-  35 

10 

304 

3b4 

-  107 

673 

34 

660 

-  26 

556 

31 

655 

404 

300 

318 

327 

346 

287 

.  288 

265 

-  284 

504 

62 

-  51 

604 

292 

335 

-  288 

333 

300 

353 

-  302 
346 

to4 

189 

-  192 

006 

-  76 

-  34 

.  44 

-  31 

106 

286 

-  300 

300 

-  289 
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hkl 

F„ 

F„ 

0 

c 

0 

0 

To6 

452 

443 

353 

424 

206 

137 

147 

177 

156 

'Z06 

228 

-  234 

230 

-  213 

306 

131 

104 

106 

114 

^6 

416 

.  406 

424 

-  392 

$06 

293 

309 

230 

277 

io6 

-  81 

76 

to6 

316 

-  304 

008 

61 

-  98 

108 

95 

147 

T08 

201 

-  236 

^8 

145 

156 

2o8 

-  80 

83 

toa 

263 

.  260 

508 

-  78 

-  4 

608 

118 

114 

708 

88 

.  126 

110 

635 

-  626 

521 

-  631 

210 

448 

429 

433 

442 

310 

201 

195 

185 

206 

410 

537 

-  510 

583 

-  551 

510 

-  115 

-  56 

62 

-  71 

610 

213 

210 

265 

239 

710 

106 

-  96 

011 

70 

-  53 

79 

-  59 

879 

-  863 

803 

-  870 

Til 

431 

438 

388 

433 

211 

169 

202 

168 

199 

111 

695 

-  716 

653 

-  713 

ill 

282 

286 

274 

296 

311 

264 

-  262 

274 

-  266 

411 

430 

-  441 

441 

-  466 

5l1 

^51 

446 

450 

458 

ill 

127 

-  144 

150 

153 

511 

139 

-  136 

141 

-  154 

611 

180 

183 

238 

212 

111 

316 

.  287 

327 

-  318 

Zii 

-  44 

-  47 

711 

150 

133 

012 

376 

-  360 

380 

-  358 

112 

619 

611 

547 

607 

112 

276 

.  302 

309 

-  297 

212 

-  80 

80 

71 

62 

I12 

853 

966 

768 

976 

312 

242 

-  243 

284 

.  269 

312 

275 

305 

283 

322 

412 

201 

204 

256 

224 

I12 

425 

.  429 

433 

.  438 

I 
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'0 

80 

79 

95 

180 

221 

194 

.  144 

185 

-  176 

294 

344 

327 

71 

.  80 

-  358 

309 

-  350 

349 

318 

346 

-  273 

283 

.  268 

.  121 

115 

.  124 

485 

380 

459 

-  394 

424 

.  412 

-  225 

203 

-  238 

124 

141 

133 

-  498 

477 

-  499 

150 

133 

283 

291 

299 

106 

-  153 

199 

230 

209 

-  214 

221 

-  243 

625 

539 

611 

-  196 

185 

-  195 

387 

397 

385 

-  94 

97 

-  94 

-  406 

441 

-  383 

280 

300 

288 

193 

177 

196 

-  35 

17 

492 

450 

501 

71 

-  80 

-  93 

62 

-  72 

-  32 

-  35 

-  17 

194 

247 

218 

265 

212 

242 

-  181 

166 

-  183 

-  56 

44 

.  62 

-  104 

106 

-  105 

-  378 

397 

-  349 

166 

185 

174 

303 

238 

298 

-  35 

-  41 

272 

238 

258 

-  332 

309 

-  332 

-  61 

219 

71 

-  58 

-  190 

168 

-  175 

183 

194 

178 

90 

79 

91 

199 

212 

195 

60. 


hkl 

F 

F 

0 

c 

0 

c 
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423 

364 

309 

351 

216 

88 

-  79 

316 

-  96 

-  79 

71 

-  50 

¥16 

-  65 

-  5^ 

-  44 

.  28 

2l6 

266 
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256 

265 

616 
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6 

716 

74 
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017 

148 
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I17 

106 

127 

217 
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417 
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.  46 

|l7 

271 

258 
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-  25 

717 
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018 

120 

135 

T18 
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-  43 

218 

.  100 

-  52 

318 

173 

176 
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164 

144 

020 
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120 
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212 

.  200 

220 
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300 
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320 
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97 
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79 
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71 
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-  437 
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409 
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462 
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621 

.  62 

21 

?21 
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221 
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721 
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-  309 

388 
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T22 

104 

89 

124 

95 
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322 
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422 

i^22 
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?22 
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023 
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123 

223 

223 
123 

323 
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423 

123 
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024 
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124 
^4 
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524 
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025 
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125 
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lf25 
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P 

P^ 

0 

c 

0 

c 

35 
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92 
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97 
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-  13 

18? 
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203 

.  214 

79 
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212 
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71 
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.  44 

22 
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-  1 

443 
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238 
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437 
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6 
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583 
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232 
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264 

300 
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197 

.  167 

221 
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556 

-  530 

539 
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230 
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212 

197 
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35 

55 

233 
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62 

70 
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79 
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70 

88 

73 
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98 
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53 

62 
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97 

86 
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-  351 

370 
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62 


hkl 


525 


026 
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730 


031 

ni 

T31 

^1 

231 

231 

331 
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F„ 

F 

F 

0 

c 

0 

c 

0120 

-  143 

.  Ill 

115 

-  113 

1120 

-  143 

98 

115 

137 

2120 

184 

154 

177 

199 

3120 

88 

-  91 

0121 

35 

-  50 

T121 

35 

69 

?121 

-  44 

26 

0122 

177 

162 

168 

204 

1122 

44 

77 

T122 

-  130 

-  78 

53 

.  81 

2122 

53 

-  94 

5122 

-  132 

-  31 

-  35 

.  22 

3122 

169 

164 

150 

212 

0123 

-  26 

38 

T123 

35 

-  97 

1131 

175 

147 

T131 

183 

-  153 
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